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ABSTRACT 



Aims. We aimed at finding a physical explanation for the occurrence of macroturbulence in the atmospheres of hot massive stars, a 
phenomenon found in observations since more than a decade but yet unexplained. 

Methods. We computed time series of line profiles for evolved massive stars broadened by rotation and by hundreds of low-amplitude 
nonradial gravity-mode pulsations which are predicted to be excited for evolved massive stars. 

Results. In general, line profiles based on macrotubulent broadening can mimic those subject to pulsational broadening. In several 
cases, though, good fits require macroturbulent velocities that pass the speed of sound for realistic pulsation amplitudes. Moreover, 
we find that the rotation velocity can be seriously underestimated by using a simple parameter description for macroturbulence rather 
than an appropriate pulsational model description to fit the line profiles. 

Conclusions. We conclude that macroturbulence is a likely signature of the collective effect of pulsations. We provide line diagnostics 
and their typical values to decide whether or not pulsational broadening is present in observed line profiles, as well as a procedure to 
avoid an inaccurate estimation of the rotation velocity. 
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1. Introduction of the phenomenon of 
macroturbulence 

Stars are gaseous bodies that transfer hydrogen into helium 
through nuclear burning in their core during the largest part of 
their lives. A variety of evolved stars results after the exhaus- 
tion of the core hydrogen burning. It is the birth mass of the star 
that determines which evolutionary path the evolved star will 
follow. Here, we are concerned with stars whose birth mass is 
above ten solar masses. Such massive stars perform subsequent 
nuclear burning cycles until their core is composed of iron, after 
which they collapse as supernova. While this broad picture of 
stellar evolution is well understood and in agreement with vari- 
ous types of observations, we still lack knowledge of important 
aspects of the physics and dynamics inside massive stars and of 
their consequences for the stellar life. 

One particular shortcoming in the description of the physics 
of stellar atmospheres of massive stars is the need to introduce 
an ad-hoc velocity field, termed macroturbulence, at the stellar 
surface in order to bring the observed shape of spectral lines 
into agreement with observations. While evidence for the oc- 
currence of such macroturbulence in hot stars was established 
since more than a decade (Howarth et al. 1997), there is still 
no physical explanation for this phenomenon available. This un- 
satisfactory situation has become ever more problematic as the 
data improved in quality in terms of resolving power and signal- 
to-noise (S/N) ratio and in quantity in terms of the number of 
stars that have been studied with high-resolution spectroscopy. 
It turns out that the macroturbulent velocities required to explain 



high-quality observations are supersonic in many of the studied 
stars, which would point to highly dynamical atmospheric mo- 
tion whose cause is unknown (Ryans et al. 2002, Lefever et al. 
2007, Markova & Puis 2008). Here, we provide a natural phys- 
ical explanation for this phenomenon in terms of the collective 
effect of numerous stellar pulsations of low amplitude. 

Velocity fields of very different scales occur in the atmo- 
spheres of stars. Apart from the rotational velocity which can 
vary from zero speed up to the critical value, line synthesis codes 
also include a certain amount of microturbulence (of order a few 
kms"') to bring the observed profiles in the spectra of stars into 
agreement with the data. Microturbulence is defined as a phe- 
nomenon related to velocity fields with scales shorter than the 
mean free path of the photons in the atmosphere (e.g.. Gray 
2005 for a thorough explanation). Microturbulence and rotation 
are usually treated as time-independent processes leading to line 
profile broadening. 

In contrast to microturbulence, macroturbulence refers to ve- 
locity fields with a scale larger than the mean free path of the 
photons (with mesoturbulence as the intermediate situation - 
e.g.. Gray 1978). Macroturbulence was mainly introduced and 
studied in the context of cool stars (e.g. Gray 1973, 1975, 1978). 
Various descriptions have been proposed in the literature (see 
Gray 2005 for an overview), among which an isotropic model 
and a radial-tangential model are the most common ones. Both 
these models will be considered here. 

Values for the micro- and macroturbulence are usually de- 
rived from line-profile fits of single snapshot spectra. Here we 
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are focused on such applications to massive hot stars, whose 
microturbulent velocities are usually below 15kms"' (e.g., 
McErlean et al. 1998; Villamariz & Herrero 2000). The pub- 
lished values of macroturbulence, on the other hand, are usually 
well above this value, reaching up to 90kms ' (Lefever et al. 
2007; Markova & Puis 2008). An important omission so far in 
the derivation of macroturbulence is that time-dependent veloc- 
ity phenomena also occur, besides rigid surface rotation and tur- 
bulence. The best known example of such a phenomenon is stel- 
lar pulsation, which causes asymmetric line-profile variations 
(e.g., Aerts & De Cat 2003). A natural step is thus to inves- 
tigate whether the needed macroturbulence may be cormected 
with the omission of pulsational broadening in the line synthe- 
sis codes used for fundamental parameter estimation. In fact, for 
pulsating stars along the main sequence, one also needs to add 
some level of macroturbulence whenever one ignores (some of) 
the detected pulsations in line-profile fitting of time-resolved or 
averaged spectra (e.g.. Morel et al. 2006). We investigate this 
hypothesis in the present paper. 

2. Computations of pulsationally broadened 
spectral line profiles 

Massive stars are exposed to pulsations during several phases 
of their life. On or near the main sequence, these pulsations are 
usually driven by a heat mechanism acting in the metal opac- 
ity bump at a temperature near some 200,000 degrees (e.g., Cox 
et al. 1992, Pamyatnyh 1999). In the recent and rapidly grow- 
ing research field of asteroseismology, observed pulsations are 
exploited by scientists to probe the badly known physical pro- 
cesses inside stars (e.g., Cunha et al. 2007, Aerts et al. 2009), 
just as it was done in helioseismology for the Sun (e.g., Gough 
et al. 1996). Asteroseismology was proven to be a valid tool to 
study the interior of massive main-sequence stars (e.g., Aerts 
et al. 2003) and may imply a unique opportunity to probe the 
internal layers, including the deep convection zone around the 
hydrogen-burning shell, of evolved stars as well. The discovery 
of gravity-mode pulsations in the Bllb supergiant HD 163899 
from spacebased high-precision photometry measured with the 
Canadian space mission MOST (Saio et al. 2006) and in a sam- 
ple of 40 B supergiants (Lefever et al. 2007) are steps in this di- 
rection. We refer the reader to Aerts et al. (2009) for a thorough 
description of stellar pulsation in all of its aspects, including the 
particular properties of the eigenfrequencies and eigenfunctions 
of pressure and gravity modes. 

It is well known that stellar pulsations imply a time- 
dependent variation of the shape of spectral lines (e.g., Aerts 
& De Cat 2003 for a review. Chapters 4 and 6 of Aerts et al. 
2009). Despite this, the estimation of the rotational and macro- 
turbulent velocity in evolved massive stars have so far usually 
been done from a single snapshot of the stellar spectrum, and 
assuming that no time-dependent phenomena are present. Here, 
we investigate to what extent stellar pulsations affect the estima- 
tion of the surface rotation and macroturbulent velocities when 
ignoring pulsations, as was often done so far in the literature. 
Hereto, we computed numerous sets of line profiles due to pul- 
sations expected in B supergiants. 

2.1. Input for the simulations 

Simulating line profile variations due to excited oscillations of a 

star requires the following steps: 

1. the computation of an equilibrium stellar structure model; 



2. the computation of the excited oscillation frequencies of the 
stellar model; 

3. the computation of the oscillation eigenfunctions in the line- 
forming region of the stellar atmosphere; 

4. the computation of the observed line profile as seen by a dis- 
tant observer, whose line of sight is inclined with respect to 
the symmetry axis of the oscillations. 

For points 1 and 2 we considered a realistic case and com- 
puted a stellar evolution model which is representative for the 
evolved pulsating Bllb star HD 163899 with the Code Liegeois 
d'Evolution Stellaire (Scuflaire et al. 2008). This model has 
the following parameters: T^s = 18, 200 K, log g = 3.05, 
RIRq = 17.8, XogiLILo) = 4.5, M/Mq = 13, Z = 0.02 and 
an age of thirtheen milhon years. It approximates well the po- 
sition of HD 163889 in the Hertzsprung-Russell diagram (Saio 
et al. 2006). We determined its excited pulsation modes of az- 
imuthal order zero with a non-adiabatic pulsation code MAD 
(Dupret 2001). We considered all modes up to degree ten, as 
it is well-known that partial geometrical cancellation effects in- 
crease with increasing mode degree (e.g.. Chapter 6 of Aerts et 
al. 2009 for a full description of these eff'ects) and, moreover, 
we needed to keep the computation time feasible. We found 241 
modes with degree { from 1 to 10 to be excited. All of them 
are gravity modes, with frequencies ranging from 0.08 to 0.68 
cycles per day, a result typical for gravity modes. 

Regarding points 3 and 4, it was shown by De Ridder et al. 
(2002) that the temperature and gravity variations in the Une- 
forming region due to the pulsations do not affect the line pro- 
file variations of a non-rotating star appreciably. This conclusion 
was based on the computation of temperature and gravity vari- 
ations for the stellar interior and for the atmosphere, and apply- 
ing a matching in a connecting layer which separates the region 
where the diffusion approximation breaks down from the stellar 
interior where it is valid (Dupret et al. 2002). This justifies the 
use of the basic line profile theory as described in Aerts et al. 
(1992). That framework allows the computation of the velocity 
eigenvectors of the modes in the line-of-sight for a linear Umb- 
darkening law. It makes use of the velocity perturbations in a 
single line-forming layer of the atmosphere to predict the line 
profile variations, while ignoring temperature, gravity, and rota- 
tional effects. Ideally, rotational effects should be included in the 
computations, given that the ratio of the rotation to pulsation fre- 
quencies can be of order one. Theories including a non-adiabatic 
treatment of rotational effects due to the Coriolis force are avail- 
able for the stellar interior, where the diffusion approximation 
is valid (e.g., Lee 2001, Townsend 2005). A study as the one 
by Dupret et al. (2002) which treats the velocity and tempera- 
ture perturbations of a rotating star in the very outer atmosphere 
is not yet available. Developing it is beyond the scope of the 
present work, which is simply to generate profiles due to pul- 
sations with properties similar to those observed and interpret 
them as macroturbulence. 

Rotation splits the frequencies of modes into 2£+\ multiplet 
components (e.g., Aerts et al. 2009). As there is currently no 
theory to provide us with the excitation of rotationally spUtted 
modes, nor with the amplitudes of the modes, we assume that all 
modes with azimuthal orders m ranging from -{ to t are excited 
with equal amplitudes in the line-forming region and we assume 
these amplitudes to be oc {( + 1) ' in the notation by Aerts 
et al. (1992). We checked that changing these assumptions does 
not alter the conclusions presented here, by considering also the 
case where only axisymmetric or sectoral modes would be ex- 
cited and by choosing different amplitude laws. Provided that a 



Aerts, Puis, Godart, Dupret: Macroturbulence in hot massive stars 3 




-100 100 -100 100 



V (km s ^) V (km s ^) 

Fig. 1. Noiseless pulsationally and rotationally broadened profiles (thin lines) are compared with the profile without pulsational 
but with rotational broadening (dashed line). The input parameters are the stellar inclination angle, the amplitude of the individual 
modes, and the projected rotation velocity, (/, a, v sin/). Their values are as follows: a: (60°, 1.0, 25), b: (60°, 1.0,45), c: (20°, 1.0,45), 
d: (60°, 0.5, 65), e: (60°, 0.5, 85), f: (60°, 0.2, 125), where a and v sin / are expressed in km s"' (see text for a definition of a). Note 
that the line features in panels d, e, f become less visible as the S/N ratio decreases. They essentially disappear for S/N < 100 (see 
also Fig.[3j. 



sufficient number of modes are included in the line broadening 
computations (typically at least a few hundred), our conclusions 
remain the same and are thus independent of the adopted am- 
plitude distribution. The conclusions are not dependent on the 
particular stellar model either 

In total, the 241 excited m - Q modes give rise to 2965 
multiplet components. We computed the collective effect of all 
these 2965 gravity modes on simulated line profiles. The simula- 
tions were made such as to mimic the effect on the Si III 4553 A 
line in the spectrum of a star with the fundamental parame- 
ters of HD 163899. In our simulations, we approximated the 
SiIII4553A Une by a Gaussian profile of width lOkms ' and 
equivalent width of 0.25 A, and we adopted a linear limb dark- 
ening law with a fixed coefficient equal to 0.364 across the line. 
These values were also fixed when computing the fits to the pul- 
sationally broadened profiles. In this way, we are sure that our 
conclusions on the macroturbulence are not affected by adopt- 
ing a wrong microturbulence or by a varying limb darkening 
coefficient across the spectral line. We limited ourselves to tun- 
ing towards this one Si spectral line, since (i) it is an important 
diagnostical line of intermediate strength, (ii) it is (almost) not 



contaminated by non-Gaussian broadening (such as Stark broad- 
ening in the case of hydrogen and helium lines) or wind effects, 
and (iii) it is the line selected for almost all of the pulsating early 
B stars so far as it turned out to be best suited to derive their pul- 
sation characteristics (Aerts & De Cat 2003). On the other hand, 
our approximation of a constant Gaussian intrinsic line implies 
that our analysis is valid for any metal line of this width in the 
spectrum. We computed time series of profiles for 50 timings 
taken from a concrete line profile study (Stefl et al. 1999) with a 
total time span of 65 days. 

We simulated various time series of 50 profiles each, tak- 
ing into account pulsational and rotational broadening, besides 
the intrinsic broadening of the spectral line. In our computations, 
we considered five values of the projected rotation velocity v sin / 
(25,45,65, 85, 125kms"'). We limited to this range of vsin/, for 
which the equatorial rotation velocities remain below 50% of the 
critical velocity of the stellar model (305 km s"' ). As can be seen 
in Fig. 1 of Aerts et al. (2004), this implies sufliciently small rel- 
ative changes of the local radius, temperature, gravity, and lumi- 
nosity to ignore the centrifugal force in the computation of the 
equilibrium structure model of the star 
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Fig. 2. Distribution of the projected pulsational velocity over the 
stellar surface as measured by a distant observer whose line-of- 
sight is inclined by 60° with respect to the rotation axis of the 
star, for four distributions of the pulsational amplitudes Vp (see 
text for further explanation). 



We adopted an inclination angle / between the rotational 
axis and the line-of-sight of 60°, but we also considered 20° 
for the case of vsin/ = 45kms The symmetry axis of the 
pulsations was taken equal to the rotational axis, as is usually 
done in pulsation studies of non-magnetic stars. Regarding the 
pulsational broadening, we considered four distributions for the 
intrinsic amplitude of the modes in the line-forming region: 
Vp = al{l + 1) with a - 1.0,0.5,0.2,0.1kms"', again using 
the notation of Aerts et al. (1992). This means that the radial 
component of the pulsational velocity vector is proportional to 
Vp while the transversal component is proportional to v^K with 
K = GM/aP'R^ with G the gravitational constant, M and R the 
mass and radius of the star, and o) the angular frequency of the 
mode (see, e.g., Aerts et al. 2009). For the adopted model we 
consider here, the ^T-values of the considered modes range from 
0.3 to 25. The choice of these amplitude distributions was made 
to end up with a realistic peak-to-peak variation of the radial 
velocity as in published observed time series of the few super- 
giant B stars for which such data are available - see Figs 5 and 
6 in Kaufer et al. (1997), Figs 2 and 3 in Prinja et al. (2004) and 
Fig. 2 in Markova et al. (2008). These studies have led to radial- 
velocity variations with peak-to-peak amplitudes between 5 and 
20 km s"' . Our amplitude distributions for Vp were taken accord- 
ingly, i.e., the collective effect of the 2965 gravity modes with 
the amplitude distributions we adopted results in radial-velocity 
variations similar to the observed ranges (see (v) in Fig. |6] dis- 
cussed in Sect. |2.2t . In this way, we are sure to have generated 
realistic profile variations, irrespective of the limitations of the 
line profile theory discussed above. A summary of the input pa- 
rameters of the simulated line profile sets, along with some of 
their computed quantities discussed below, is given in Table[T] 

The radial velocity is an integrated quantity over the stellar 
surface. Pulsating stars have time-dependent asymmetric line- 
profile variations. It is common to characterise the line profile 
shapes by their three lowest-order moments, which represent the 
centroid velocity (v), the width (v^) and the skewness (v-*). A 
practical guide to compute these quantities, as well as their for- 
mal definition in terms of the surface velocity eigenfunctions, 
is provided in De Ridder et al. (2002) and more extensively in 



Chapter 6 of Aerts et al. (2009). Aerts et al. (1992) and Aerts 
(1996) provided a thorough discussion of these quantities and 
their suitability to interprete them in terms of pulsation the- 
ory, thus allowing an identification of the spherical wavenum- 
bers ({, m) from observed time series of moment variations. We 
computed these three quantities for the simulated line profiles, 
mainly to show their relation with the derived macroturbulent 
velocity values obtained when ignoring the pulsational broaden- 
ing, as will be discussed in Sect. l2.2l The first moment (v) is the 
radial velocity of the star, integrated over the stellar disc, with re- 
spect to the centre of mass of the star (i.e. it varies around value 
zero during the pulsation cycle); it is thus directly comparable to 
the measured radial-velocity variations reported in the literature 
which are usually based on Gaussian fits to the profiles. 

In order to end up with peak-to-peak radial-velocity vari- 
ations of order 20 km s"\ as measured for several supergiant 
B stars from metal lines, numerous of the individual surface 
elements must experience a far larger individual pulsation ve- 
locity. In the case of radial pulsations and in the approxima- 
tion of the adopted linear limb darkening law, this means that 
the entire surface moves up and down with a velocity of about 
20 X 1.5 = 30km s"'. This also implies that measured radial- 
velocity variations above typically 40 km s"' are the results of 
shock phenomena in the atmosphere of radial B-type pulsators, 
leading to a so-called "stillstand" in the radial-velocity curve in 
the case of radial modes. Examples can be found in Aerts et al. 
(1995), Saesen et al. (2005) and Briquet et al. (2009) for the 
^Cep stars BW Vul, ^' CMa and V1449Aql, respectively. Such 
stillstand was so far not observed for B supergiants, so we ex- 
pect the majority of the surface elements to move subsonically 
(which does not imply that some elements may encounter super- 
sonic speeds). 

For non-radial gravity modes, a wide variety of surface ve- 
locities occurs across the stellar surface and shock phenomena 
are much harder to detect in integrated quantities, such as mo- 
ments or equivalent widths. We show in Fig. |2] the distribution 
of the line-of-sight components of the total pulsational velocity 
vectors, which result from the addition of all the individual vec- 
tors of the 2965 modes, for each of the surface elements (denoted 
as |vpuis(/?, 0, (p, t)\), for the four amplitude sets corresponding to 
the four a-values. The addition of the numerous eigenvectors can 
result in positive or negative projected velocity, depending on the 
phases of the modes and on the place (R, 6, <p) on the surface. We 
expect that, in most of the surface points and for most of the tim- 
ings, positive and negative contributions tend to lead to a limited 
value of the overall pulsation velocity due to cancelling of posi- 
tive and negative mode velocities, since we assumed there to be 
no phase relation between the modes. It can be seen from Fig.|2] 
that, for all four amplitude sets, the majority of the surface ele- 
ments indeed are seen to move subsonically. For a = l.Okms"', 
supersonic speeds in the line-of-sight are encountered for a con- 
siderable fraction of the surface elements, but still in less than 
half of them such that the radial-velocity variations remain be- 
low 20km s"' (see Fig. |6] discussed in Sect. 12. 2b . The adopted 
amplitude distributions thus lead to realistic peak-to-peak am- 
plitudes for the radial velocity. In this way, we are sure not to 
overestimate the effects of pulsations on the derivation of the 
macroturbulent velocity values. 

We considered profile sets without noise and with white 
noise resulting in S/N ratios of 200 and 500. This brings the 
total number of simulates profiles to 3600 (50 timings, 6 com- 
binations (v sin /, /), 4 amplitude distributions and 3 noise lev- 
els). Examples are provided in Fig.[T]and lead to the conclusion 
that some of the simulated profiles are considerably affected by 
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the collective effect of the gravity modes. In particular, the line 
wings are broader than those that would occur for a star that does 
not have pulsations. 

2.2. Estimation of the macroturbuient veiocities 

Various possibilities to describe macroturbulence have been pre- 
sented in the literature. We refer to Gray (2005) for a thorough 
discussion. In this work, we considered an isotropic macroturbu- 
lence described by a Gaussian velocity distribution (denoted as 
Aiso), as well as an anisotropic description for which the radial 
and tangential velocity fields in general have a different ampli- 
tude denoted as Ar and Aj (a so-called radial-tangential model 
— see Eq. (17.6), p.433 in Gray 2005). For the anisotropic model 
fits, we considered the two extreme cases of allowing Ar to be 
free while Ax = and Ar = while Aj was allowed to take any 
value. In this way, each of the three models for the macroturbu- 
lence is described by one free parameter 

For all the simulated profiles, we determined v sin / and the 
macroturbulence Vmacro, while ignoring the presence of pulsa- 
tional broadening, as is done in the literature, by adopting a 
goodness-of-fit approach. The normalized profiles broadened 
by both rotation and gravity-mode pulsations are denoted by 
(Aj,pi(Aj)) and those broadened by rotation and macroturbu- 
lence by (Aj, p2iAj)), with j - 1, . . .,N an index labelling the 
velocity pixels within the profile. For the computation of p2 we 
considered each of the three options Aiso, Ar, and At. Each of 
the profiles pi and p2 were given the same equivalent width. We 
computed the line deviation parameter, X, based on the classical 
statistical technique of residuals: 



2(v sin /, Vmacio) = 



A' 



1 2 



(1) 



This quantity is the standard deviation of the residual profile 
\p\ - Pil, averaged over all velocity pixels in the line profile, 
expressed in continuum units. It is thus a measure of the fit qual- 
ity, directly interpretable in terms of the S/N ratio of measure- 
ments. The optimal choice of the parameters (v sin /, vmacro) is 
then found by carefully screening the 2-dimensional parameter 
space in v sin ; and Vmacro (in steps of 1 km s"' for each of v sin / 
and Vmacro) and by identifying Em s min(,,sin/,i„,,,„o) where 
Vmacro Can be any of Aiso, Ar or Aj. Moreover, we allowed two 
options to fit the pulsationally broadened profile p i : the one such 
that the wavelength position of the minimum of p i and p2 coin- 
cide and the one such that their first moments (v) are in best 
agreement. Note that these two options are equivalent only in 
the case of symmetric profiles. For each generated profile, we 
performed these six fitting exercises, keeping in each case the fit 
with the lowest Zm(v sin /, Vmacro) as the best one. We then com- 
pared the parameters (vsin/, Vmacro) of the best fitting profile p2 
with the input vsin / of pi and analysed the values of Vmacro- Six 
prototypical examples of best fits are shown in Fig. [3] 

In Fig.|4]we show the outcome of the fit to the 3600 simu- 
lated profiles, for the three models we considered for the macro- 
turbulence. For all three models, it was found that the inclusion 
of an ad-hoc macroturbulence parameter leads to better fits than 
those obtained when only allowing rotational broadening, which 
has to be the case given that there is one more degree of free- 
dom. This is visible from Fig.|5]where we show the distribution 
of Em deduced from fits with and without allowing a parame- 
ter for macroturbulence. Figure|5] contains all simulated profiles; 
these two global distributions are the same as those for the five 
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Fig. 5. Distribution of Em, which is a measure of the fit quality 
averaged over the line profile and expressed in continuum units - 
see Eq. ([TJ, for 3600 fits without macroturbulence (full line) and 
with macroturbulence (either AjscAr, or At whichever gave the 
lowest Em, dotted line). 



separate values of v sin /, which is as expected given that the pul- 
sational broadening was simply added to the rotational broad- 
ening without any coupling between the two. One would typi- 
cally improve the fit quality from eye inspection by incorporat- 
ing macroturbulence for Em > 0.008; this corresponds with the 
dotted lines in panels b,d,e,f in Fig. [3] whose counterpart with 
macroturbulence represented by the dashed lines imply a notica- 
ble reduction in Em (for comparison, the dotted lines in panels 
a and c have Em = 0.0058 and 0.0056, respectively, and would 
probably not give rise to the introduction of macroturbulence). 
It can be seen from the distribution of Em in Fig. |5] that the fit 
quality is very good for the large majority of profiles when al- 
lowing for macroturbuient broadening. In 89% of the cases, the 
fit with macroturbulence has Em < 0.008. If we do not include 
macroturbulence, 59% of the fits has Em < 0.008. 

Returning back to Fig.|4] we deduce that the lowest Em 
values were reached for Aiso,Ar, and At in 1224, 1343, and 
1033 of the cases, so these descriptions are basically equivalent 
in appropriateness to mimic pulsational broadening. It can be 
seen from Fig. |4] that the radial model Ar needs higher values 
to achieve a good fit compared with the isotropic and tangen- 
tial model. This is logical, because the pulsation velocities of 
the modes are dominantly horizontal in nature. There are some 
differences between the fit quality in a global sense for three 
considered models, but the main conclusion is that the missing 
broadening caused by the pulsations can often only be compen- 
sated by quite large values of the macroturbulence. 

In Fig.|6]we show the value of the macroturbulence with the 
lowest Em, along with the first and third velocity moment, of the 
simulated profiles. As Fig.|4] this plot illustrates that the value 
of the macroturbulence can be very large, compatible with what 
is found in the literature, if one ignores pulsational broadening, 
even though the centroid velocity variations (v) induced by the 
pulsations are modest. The reason for this is that the line width 
is a function of the square of the velocity, and the line skewness 
is represented by (v^). Thus, one needs to compensate the line 
width and line wing shape by a large value for the macroturbu- 
lence whenever one wants to achieve a good profile fit. 
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The reported absence of line asymmetries in the literature 
must be compatible with our physical line-broadening model. 
First of all, most spectroscopic studies in which symmetric pro- 
files are mentioned rely on visual inspection of only one spec- 
trum, while line profile variability (and thus asymmetry) is al- 
most always found when multiple-epoch observations are taken. 
Visual inspection of the fits in the top and bottom of the right 
column in Fig. [3] reveals line asymmetry from one snapshot 
spectrum, while the other four profiles might give the impres- 
sion of being symmetric. Typically, the profiles simulated with 
a = 1.0 km s"' would be detectable by visual inspection of the 
profiles. However, when one computes diagnostic line quanti- 
ties, it often becomes obvious that the lines deviate from sym- 
metry even if seemingly symmetric by eye. The best diagnostic 
parameters to characterise line asymmetry in the case of pul- 
sations are the line moments. While line bisectors and velocity 
spans are often used in the cool star and exoplanet communities, 
such parameters are not suited to be interpreted in terms of pul- 
sational parameters while moments are (e.g., Aerts et al. 1992, 
Dall et al. 2006, Hekker et al. 2006). In practice one can use the 
property that the odd moments of a symmetrical line profile are 



zero. Thus, the values of (v) and (v-') of metal lines measured 
with a high resolution and high S/N ratio are well suited to de- 
cide if an observed line profile is subject to time-dependent line 
asymmetry whenever this is not obvious from visual inspection. 
As a guide, we provide the ranges of the values of the moments 
of the generated profile sets in Table[T] The values of (v) for the 
six profiles shown as full lines in Fig. [3] are a: -1.1, b: 0.8, c: 
-0.3, d: -1.0, e: -1.4, andf: 1 .9 km s ' . The corresponding val- 
ues of (y^) are a: -1805, b: 226, c: 9964, d: -45658, e: -36972, 
and f: 12470 km^ s^^*. All these values would be zero in the case 
of symmetrical profiles subject to white noise. The deviation is 
small, of order a few kms"', for (v), because this quantity mea- 
sures the centroid of the line and thus is independent of v sin / 
and the microturbulence, while these two quantities do affect 
(v^) (see Aerts et al. 1992). 

The use of the odd moments for asymmetry detection has to 
be done from a few spectra spread at least a few days in time, be- 
cause line blending of course also causes a deviation from sym- 
metry. Such deviation is time-independent, though, while the 
signature of pulsational broadening is always time-dependent 
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and has typical periodicities of several hours to a few days in 
hot massive stars. 



2.3. Consequences for the rotational velocity estimate 

An important conclusion based on Fig.|6]is that the inclusion of 
macroturbulence to obtain a line fit may result in a serious un- 
derestimation of the true projected rotational velocity, irrespec- 
tive of which description for Vmacro is used. We encountered mis- 
matches compared with the input vsin/ above 100 km s"^ The 
question thus arrises if it is not wiser to exclude a vmacro parame- 
ter to search the best value of v sin / or to resort to other methods 
to determine this parameter 

The rotation velocities derived from a fit to the pulsationally 
broadened profiles with the inclusion of macroturbulence (we 
selected the version of Aiso, Ar, and Aj which led to the low- 
est for the plot) and without it are compared with the input 
V sin / in the left panel of Fig.|2l It can be seen that the mismatch 
of the rotation velocity ranges from -20 to 40 km s" ' for fits with- 
out allowing a parameter for macroturbulent broadening while it 
ranges from -120 to 40 km s"' if a parameter for macroturbu- 
lence is allowed. From this we conclude it is better to avoid the 
inclusion of macroturbulence in a goodness-of-fit approach as 
in Eq. (Q when the goal is to achieve a good estimate of v sin /. 



Even in that case, v sin / may be quite wrong when derived from 
profiles which are pulsationally broadened. 

In view of the importance of having an appropriate v sin / es- 
timate, we also resorted to the popular Fourier Transform (FT) 
method. This method was introduced by Gray (1973, 1975). It 
was evaluated specifically for hot massive stars by Simon-Diaz 
& Herrero (2007). It allows to estimate v sin / from the first mini- 
mum of the FT of a line profile. What is often forgotten, however, 
is that its basic assumption is that the line profiles are symmetric, 
which is not the case when pulsations (or other phenomena like 
spots) occur (e.g.. Smith & Gray 1976). We thus investigated 
how robust the method is when this condition is not met. We ap- 
plied the method to all the 3600 pulsationally broadened profiles 
and derived v sin / by careful visual inspection of their FTs. A 
few of the FTs are shown in Fig. |9] while the global mismatch 
in V sin / is compared with the one obtained from the goodness- 
of-fit method for vmacro = in the right panel of Fig.|2] We 
see that the FT method outperformes the goodness-of-fit method 
when v'niacio is allowed for On the other hand, the ability of the 
FT method to estimate v sin / is also affected by the pulsational 
broadening for a fraction of the simulated profiles and also leads 
to too low estimates for vsin/. The offsets between the input 
value of V sin / and its estimate from the goodness-of-fit method 
with Vmacro = on the one hand, and from the FT method on 
the other hand, are above 10km s"' in 16% and 26% of the 3600 
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cases, respectively. For mismatches above 20km s" these num- 
bers decrease to 6% and 12%, and above 30km s"' a further de- 
crease to 1.7% and 5% occurs. 

As illustrated in Fig.|9]for a few cases, the results of the FT 
method improve appreciably when applied to the best fit of the 
line profile including only microturbulent and rotational broad- 
ening, i.e., without allowing macroturbulence. The reason is that, 
in this case, one approximates the true pulsationally broadened 
asymmetric profile by one which is symmetric and has less ex- 
tended wings such that the basic assumption of the FT method is 
fulfilled. It was already emphasized long ago by Mihalas (1979) 
that the FT method has limitations of applicability when various 
broadening functions are convolved and result in skew profiles, 
which is the situation we encounter here for the gravity modes. 
The FT method is reliable in filtering out the value of v sin / from 
the observed spectral lines, when the rotational broadening is 
very dominant while the pulsational amplitudes are very low (as 
in panel c of Fig.|9]l or when broadening due to spots or pressure 
modes occurs, which leave the line wings almost unaltered and 
affect mainly the central parts of the lines. 

In Fig. [8] we compare the input value of v sin / with the value 
deduced from a fit with and without allowing a parameter for 
macroturbulence; we also show the corresponding values of 



I'macro for each of the simulated profiles. We find an overesti- 
mation of the rotation from Em for low input v sin /, because we 
need to compensate for the pulsational broadening and this can 
only be achieved by fitting a profile with a higher v sin / than the 
input value. When allowing for a macroturbulence, however, we 
cover the entire range of projected rotation velocities between 
zero and values up to some 20 km s"' above the input value of 
V sin /, i.e., for several cases a serious underestimation of the true 
rotational velocity occurs. This mismatch increases with the in- 
put V sin / and occurs whenever the wings of the profiles are sev- 
eraly broadened due to a positive interference of the modes with 
the largest horizontal velocity amplitude at some timings in the 
beat cycle and/or when large line asymmetries occur (see, e.g., 
panels d, e, f of Fig.[3]which typically have large values of Kv-')!). 
The right panel of Fig. [8] shows that large values of Vmacro occur 
for all input values of vsin/, but the most extreme values for 
v'macio occur typically for the broader profiles due to rotation and 
the larger pulsational amplitudes (see also Fig.|4|l. 

We come to the important conclusion that, in the case of line 
profile broadening due to gravity modes, v sin / estimates are best 
derived from a simple goodness-of-fit to observed profiles, in- 
cluding only microturbulence and rotational broadening and no 
macroturbulence . 
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Table 1. Ranges of the intervals for the macroturbulence and moments of the 72 sets of simulated Une profiles with pulsational 
broadening (3 values of the S/N for each of the combinations (/, v sin /, a), with Vp = ajit + 1). See text for further explanation. 
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Fig. 7. Left: the estimated minus input rotation velocity from a fit without macroturbulence as a function of a fit with macroturbu- 
lence, for the simulations described in the text. The two different symbols indicate simulations for the four amplitude distributions 
a - 0.1, 0.2, 0.5(o), 1.0(+)kms"' as explained in the text and listed in Table[T] Right: the estimated minus input rotation velocity 
from a fit without macroturbulence as a function of the value derived from the Fourier Transform method. 



3. Implications 

The idea that macroturbulence originates from stellar pulsations 
is not new. In fact, Lucy (1976) akeady suggested pulsations as a 
potential possibility to explain macroturbulence. Unfortunately, 
he did not have the observational capabilities nor the theo- 
retical development to study the effects of pulsations on line 
profiles. Recent observations of massive stars with the CoRoT 
space mission indeed reveal the occurrence of hundreds of pul- 
sation modes with white-light amplitudes in the range 0.01 and 
0. 1 mmag which went unnoticed in ground-based data (Degroote 
et al. 2009a,b). Moreover, the discovery of massive pulsators in 
low-metallicity environments (e.g., Kolaczkowski et al. 2006, 
Narwid et al. 2006, Sarro et al. 2009) also shows that current 



excitation computations (Miglio et al. 2007) still underestimate 
the number of excited modes. Our results are thus also relevant 
for evolved stars in the Magellanic Clouds. 

As an important side result of our study, we conclude that 
the rotational velocities of evolved massive stars can be seriously 
underestimated by using line profile fits based on a description 
in terms of macroturbulence. Ironically, this finding is just oppo- 
site to previous arguments that, by neglecting macroturbulence, 
the derived v sin ; values are significantly overestimated. In or- 
der to avoid erroneous estimates of v sin /, we advise to compute 
the moments of the line profiles as well as to compare the val- 
ues of V sin / from fits with and without allowing macroturbu- 
lent broadening both by a goodness-of-fit approach and by the 



10 



Aerts, Puis, Godart, Dupret: Macroturbulence in hot massive stars 




50 75 100 

vsini (input) (km s ^) 



50 75 100 

vsini (input) (l<m s ') 



Fig. 8. Left: v sin / for fits without macroturbulence (circles) and with macroturbulence (plus signs) as a function of the input v sin /, 
for all the simulated line profiles. The circles were artificially shifted to slightly lower input v sin / for visibility purposes. Right: the 
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Fourier method. In this way, the probability of a wrong v sin / 
estimate is relatively low. 

It is remarkable that the link between pulsational broadening 
and macroturbulence, and its effect on the derivation of v sin /, 
was never thoroughly investigated, particularly since the sur- 
face rotational velocity derived from line profile fitting consti- 
tutes a crucial stellar parameter which is used to evaluate stel- 
lar evolution theory. Several authors, among which Hunter et al. 
(2008), claim to have found too low observed rotational veloci- 
ties for evolved massive stars compared with theoretical predic- 
tions. Our physical model of collective pulsational broadening 
may help resolve this discrepancy. Accurate derivations of the 
rotational velocity of massive stars are also relevant in the con- 
text of Gamma-Ray-Burst progenitor studies (e.g., Yoon et al. 
2006). In view of our results, we strongly advise to use multi- 
epoch observations, because that is the best way to estimate the 
effect of pulsational broadening. One should attempt to take at 
least ten spectra with a resolution above 30,000 and a S/N ratio 
above 200, spread over different nights, and consider the broad- 
ening of different metal lines, to achieve a valid estimate of the 
surface rotation. 

Our present study was based on simulations in which we 
considered pulsational line broadening due to velocity pertur- 
bations ignoring the Coriolis force, which were then interpreted 
as macroturbulence. The resulting simulated profiles were con- 
structed in such a way as to lead to realistic radial-velocity vari- 
ations. It might be worth to investigate how the inclusion of the 
collective effect of non-adiabatic temperature and gravity varia- 
tions in the line-forming region of a star subject to the Coriolis 
force will affect the line wing broadening and its interpretation 
in terms of macroturbulence. Irrespective of the limitations of 
present line profile theory, our conclusion is clear: ignoring time- 
dependent pulsational line broadening in line profile fits of snap- 
shot spectra may lead to the need to introduce an ad-hoc veloc- 
ity field to account for the missing broadening in the line wings. 
This implies the risk of a wrong estimation of the projected ro- 
tational velocity of the star 
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